February 2, 2008 9:46 Proceedings Trim Size: 9in x 6in 



Blas'Korea 



COOLING OF NEUTRON STARS WITH COLOR 
SUPERCONDUCTING QUARK CORES * 



D. BLASCHKE 

Fachberetch Physik, Universitdt Rostock, D-18051 Rostock, Germany 
oliubov Laboratory for Theoretical Physics, JINR Dubna, 141980 Dubna, 

Russia 

email: david@thsunl.jinr.ru 

D. N. VOSKRESENSKY 

Theory Division, GSI mbH, D-64291 Darmstadt, Germany 
Moscow Institute for Physics and Engineering, 115409 Moscow, Russia 

H. GRIGORIAN 

FB Physik, Universitdt Rostock, D-18051 Rostock, Germany 
Department of Physics, Yerevan State University, Alex Manoogian Str. 1, 
375025 Yerevan, Armenia 



We show that within a recently developed nonlocal, chiral quark model the critical 
densities for a phase transition to color superconducting quark matter under neu- 
tron star conditions can be low enough that these phases occur in compact star 
configurations with masses below 1.4 Mq. We study the cooling of these objects 
in isolation for different values of the gravitational mass and thus different compo- 
sition and structure of the interior. Our equation of state allows for a 2SC phase 
with a large quark gap A ~ 100 MeV for u and d quarks of two colors, a nor- 
mal quark matter phase and their coexistence in a mixed phase within the hybrid 
star interior. We argue that, if the phases with unpaired quarks were allowed, the 
corresponding hybrid stars would cool too fast to describe the neutron star cool- 
ing data existing by today. We incorporate other attractive channels permitting 
a weak pairing of the residual quarks which remained unpaired in the 2SC phase 
and demonstrate that the model does not contradict the cooling data if the weak 
pairing gaps are of the order of 0.1 MeV. 
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1. Introduction 

In the recent paper ^, hereafter BGV, we have reinvestigated the cooHng 
of neutron stars (NS) within a purely hadron model, i.e., ignoring the pos- 
sibility of quark cores in NS interiors. We have demonstrated that the NS 
cooling data available by today can be well explained within the "Nuclear 
medium cooling scenario", i.e., if one includes medium effects in the emis- 
sivity and takes into account a suppression of the 3^2 neutron gap. In 
a subsequent work ^ we have shown that this result does not exclude the 
possibility that neutron stars might possess large quark matter cores that 
extend up to more than half of the star radius. Such a hybrid structure 
gives room for a whole variety of additional scenarios of compact star cool- 
ing which fall into two classes: either nuclear and quark matter phases have 
similar cooling behavior (homgcneous cooling) or the faster cooling of the 
one phase is compensated by the slower cooling of the other (inhomoge- 
neous cooling). In the present contribution we will report on our results 
within the former, homgeneous cooling scenario of hybrid stars and what 
implications the comparison with present-day cooling data may provide for 
the EoS and transport properties of quark matter. 

2. Color superconductivity 

The quark-quark interaction in the color anti-triplet channel is attractive 

driving the pairing with a large zero-temperature pairing gap A ~ 100 MeV 
for the quark chemical potential jig ^ (300-^500) MeV, cf. ^-^j see review ^ 
and Refs therein. The attraction comes either from the one-gluon exchange, 
or from a non-perturbative 4-point interaction motivated by instantons ^ 
or from non-perturbative gluon propagators ^. Various phases are possi- 
ble. The so called 2-flavor color superconductivity (2SC) phase allows for 
unpaired quarks of one color, say blue. There may also exist a color-flavor 
locked (CFL) phase ^ for not too large values of the dymanical strange 
quark mass or other words for very large values of the baryon chemical 
potential ^, where the color superconductivity (CSC) is complete in the 
sense that the diquark condensation produces a gap for quarks of all three 
colors and flavors. The value of the gap is of the same order of magnitude 
as that in the two-flavor case. There exist other attractive quark pairing 
channels for quarks that can't participate in 2SC and CFL pairing. These 
weak pairing channels are characterized by gaps typically in the interval 
~ 10 keV -^1 MeV, and a prominent example which will be used in the 
present work is the spin-1 pairing channel of single color diquarks in the 
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isospin singlet state, also named as color-spin locking (CSL) phase. 

The high-density phases of QCD at low temperatures may exist in the 
interiors of hybrid stars affecting their cooling, rotation and magnetic field 
evolution, cf. 11.12,13,14,15 

3. Hybrid stars 

In describing the hadronic part of the hybrid star we exploit the Argonne 

V18 + 6v + UIX* model of the EoS given in 1^, which is based on the most 
recent models for the nucleon-nucleon interaction with the inclusion of a 
parameterized three-body force and relativistic boost corrections. Actually 
we adopt here an analytic parameterization of this model by Heiselberg 
and Hjorth- Jensen 1^, hereafter HHJ. The latter uses a compressional part 
with the incompressibility K ~ 240 MeV, and a symmetry energy fitted 
to the data around nuclear saturation density that smoothly incorporates 
causality at high densities. The density dependence of the symmetry energy 
is very important since it determines the value of the threshold density 
for the DU process {n^^). The HHJ EoS fits the symmetry energy to 
the original Argonne V18 + Sv + UIX* model yielding nj?^ ~ 5.19 no 
(Mfu 1.839 Mq). 

The 2SC phase occurs at lower baryon densities than the CFL phase, 
see is>i9. For applications to compact stars the omission of the strange 
quark fiavor is justified by the fact that central chemical potentials in star 
configurations do barely reach the threshold value at which the mass gap 
for strange quarks breaks down and they appear in the system 

We will focus on the model of the quark EoS developed in The 
Gaussian, Lorentzian and NJL type cutoff formfactors were studied. The 
Lorentzian interpolates between a soft (Gaussian type, a ~ 2), and a hard 
(NJL, a > 30) depending on the value of the parameter a. We will further 
work out two possibilities of the Gaussian and the Lorentzian formfactors. 

In some density interval at the first order phase transition there may 
appear the region of the mixed phase, see Ref. disregarded finite 
size effects, such as surface tension and charge screening. Refs on the 
example of the hadron-quark mixed phase have demonstrated that finite 
size effects might play a crucial role substantially narrowing the region of 
the mixed phase or even forbidding its appearance. Therefore we omit the 
possibility of the hadron-quark mixed phase in our model where the quark 
phase arises by the Maxwell construction. For the case of the Gaussian 
formfactor the quark core appears for n > iic = 0.34fm~'^ (M > 0.852 Mq), 
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Figure 1. Mass - radius relations for compact star configurations with different EoS: 
purely hadronic star with HHJ EoS (dashed line), stable hybrid stars with HHJ - SMq' 
EoS (solid line) and with HHJ - SM^^' EoS (dash-dotted line). 



for the Lorentzian formfactor n > ric — O.SOfm^^ (M > 1.81 Mq). In the 
following, we will not further discuss the Lorentzian case which gives a 
marginal quark core in a small mass range only, see Fig. 1. 

A large difference between chemical potentials of u and d quarks forbids 
the pure 2SC phase, cf. The CFL phase is still not permitted at 

such densities. Ignoring the possibility of a weak coupling we have two 
possibilities: either the quark matter for n > Uc is in the normal phase, or 
there appears a region of the 2SC - normal quark mixed phase. Following 

Ref. considered the latter possibility disregarding finite size effects 
and has found the possibility of a wide region (for Gaussian formfactor) of 
the 2SC - normal quark mixed phase instead of a pure 2SC phase. In the 
given case arguments of '^^ are relaxed since the surface tension on the 2SC 
- normal quark boundary should be much less compared to that for the 
quark - hadron boundary. Indeed, the surface tension is proportional to 
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the difference of the energies of two phases, being very small in the given 
case, oc (A//Zg)^ <C 1. 

In Fig. 1 we present the mass-radius relation for hybrid stars with HHJ 
vs. SM EoS. Two sets of configurations given by Gaussian (solid lines) 

and Lorentzian (dash-dotted lines) formfactors are stable, sec also Rcf. 
where similar results for the nonlinear Walecka model EoS (RMF) have 
been found. 

If one switches on the possibility of the weak coupling, e.g. the CSL 
pairing channel, see 10,19,20 .^^i the quarks in the normal phase may acquire 
a corresponding spin-1 pairing gap, typically A ~ 10 keV -^10 MeV. In such 
a way all the quarks may get paired, either strongly in the 2SC channel or 
weakly in the CSL one. 

4. Cooling 

For the calculation of the cooling of the hadron part of the hybrid star 
we adopt the same model as in BGV. The main processes are medium 
modified Urea (MMU) and pair breaking and formation (PBF) processes. 
The HHJ EoS was adopted. In Fig. 2 (Fig. 20 of BGV) we show the 

cooling of different NS calculated within the hadron model of BGV. We 
use a fit-law for the relation between surface and interior temperatures, 
see BGV. Possibilities of the pion condensation and of the other so called 
exotic processes are for simplicity suppressed. Direct Urea is irrelevant in 
this model up to very large NS mass M > 1.839 Mq. ISq neutron and 
proton gaps are taken the same as in paper , whereas 3P2 neutron gap is 
suppressed by the factor 0.1, see Fig. 4 of BGV. 

For the calculation of the cooling of the quark part of the hybrid star 
we are basing on the model We include the most efficient processes: the 
quark direct Urea (QDU) processes on unpaired quarks, the quark modified 
Urea (QMU), the quark bremsstrahhuig (QB), the electron bremsstrahlung 
(EB) Following we include the emissivity of the quark pair formation 
and breaking (QPFB) processes. The specific heat incorporates the quark 
contribution, the electron contribution and the gluon-photon contribution. 
In the CSL phase the specific heat is proportional to T^. This new 
term does not significantly affect the total specific heat since the main 
contribution comes from electrons. The heat conductivity contains quark, 
electron and gluon terms. 

We arc basing on the picture presented in Fig. 2 and add the contri- 
bution of the quark core. For the Lorentzian formfactor the quark core 
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Figure 2. Cooling curves according to the nuclear medium cooling scenario, see Fig. 20 
of BGV. The labels correspond to the gravitational masses of the configurations in units 
of the solar mass. 

appears only for M > 1.81 M©, see Fig. 1. The existing cooling data are 
not affected, thereby. 

For the Gaussian formfactor the quark core occurs already for M > 
0.852 Mq according to the model see Fig. 1. Most of the relevant NS 
configurations (see Fig. 2) are then affected by the presence of the quark 
core. First we check the possibility of the 2SC+ normal quark phases. 

Fig. 3 shows the cooling curves calculated with the Gaussian ansatz. 
The variation of zero temperature gaps for the strong pairing of quarks 
within 2SC phase in the interval A ~ 20 200 MeV only slightly affects the 
results. The main cooling process is the QDU process on normal quarks. 
We see that the presence of normal quarks leads to too fast cooling. The 
data could be explained only if all the masses lie in a very narrow interval 
(0.82 < M/Mq < 0.90 in our case). The existence of only a very narrow 
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Figure 3. Cooling curves for hybrid star configurations with Gaussian quark matter core 
in the 2SC phase. The labels correspond to the gravitational masses of the configurations 
in units of the solar mass. 



mass interval in which the data can be fitted seems us unreahstic as by 
itself as from the point of view of the observation of the NS with different 
masses: M ~ 1.41 Mq and M ~ 1.25 Mq, cf. Thus the data can't be 
satisfactorily explained. 

In Fig. 4 we permit the weak pairing for all quarks which were assumed 
to be unpaired. We use A ~ 1 MeV for the corresponding quark gap. Fig. 
4 demonstrates too slow cooling. 

In Fig. 5 we again allow for the weak pairing of those quarks which 
were assumed to be unpaired in Fig. 3, but now we use A = 50 keV for 
the corresponding CSL quark gap. The data are appropriately fitted. The 
" slow cooling" data are explained by the cooling of NS either without quark 
core (M ~ 0.5 Mq, see Fig. 2) or with a quark core. Although the majority 
of experimental cooling points are covered by NS masses in a very narrow 
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Figure 4. Cooling curves for hybrid star configurations with Gaussian quark matter 
core in the 2SC+CSL phase. The weak pairing gap is 1 MeV. The labels correspond to 
the gravitational masses of the configurations in units of the solar mass. 



mass interval M — 0.86 ~ 0.87 Mq the difference witli Fig. 3 is crucial. 
We see that with the higher masses, up to 1.6 Mq one may cover the rapid 
cooling point (Vela). Besides, the variation of parameters allows to shift all 
the curves up what permits to essentially broaden the mass interval that 
would cover the data. E.g. with the so called Tsuruta law for the — Tin 
relation the corresponding mass interval is M = 0.86 1.79 Mq. With 
these remarks we showed that the 2SC+CSL+hadron scenario allows to fit 
the data as well as the purely hadronic scenario, see ^. To be more realistic 
one should further inlcude density dependences of the gaps what we intend 
to do in the forthcoming publication. 



2, 2008 9:46 Proceedings Trim Size: 9in x 6in 



Blas'Korea 



9 




log(t[yr]) 

Figure 5. Same as Fig. 4 but with a weak pairing gap of 50 keV. 



Conclusion 

Concluding, we demonstrated that the present day cooling data can be 
explained not only by a purely hadronic structure of NS interiors but also 
by a hybrid one with a complex pairing pattern, where quarks are partly 
strongly paired within the 2SC channel, and partly weakly paired within 
the CSL channel with gaps A ^ 50 keV. We conclude also that our choice 
of a density-independent weak pairing gap could be the reason why the 
mass interval for explaining slow and intermediate cooling data is very 
narrow. As it is well-known that the CSL gap should have a strong density 
dependence the fastening of the cooling by increasing the star masses 
should be partly compensated for by the corresponding increase of the weak 
pairing gap. Corresponding calculations are under way. 
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